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1. Introduction 



The precise determination of the elements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
can impose important constraints on physics beyond the Standard Model (SM). One combination of 
CKM matrix elements that plays a relevant role in this analysis is ^ , which is related to Bo — #o 
mixing. 

In particular, this combination of CKM matrix elements can be extracted from the precisely 
experimentally measured quantities AM S and AMj, which are the mass differences between the 

B° systems respectively. The relation is 



heavy and light mass eigenstates in the B s 
given by 



V tl , 



B { > and B» 



V ts 




(1.1) 



The masses Mb s and Mg d , and the corresponding mass differences are known experimentally 



=, however, an accurate and consistent lattice 



with very high precision pi]. For the ratio t 

calculation that fully incorporates vacuum polarization effects is not yet available. Our goal is to 
perform such a calculation and reduce the theoretical errors in the ratio E, to a few percent. This 
will provide us with a high precision determination of the CKM ratio in ( |1.1| ). 

The products of Bo decay constants and bag parameters in ( [ET| ) are determined by matrix 
elements between Bo and Bo of the four-fermion operators appearing in the effective hamiltonian 
that describes AB = 2 processes. The non-perturbative inputs for the calculation of Ar s and Ar^ 
(with Ar the width difference between the light and heavy mass eigenstate) are also given by 
this kind of hadronic matrix elements. For completeness, we are studying all the matrix elements 
needed to make theoretical predictions for AM S , AM C /, AT S and ATj. 



2. Simulation details and milestones in the calculation 

The four-fermion operators whose matrix element between Bo and Bo are needed to make a 
complete study of Bq and Bq mixing in the SM are 

013 ee [bW] V _ A [b j q j ] V _ A ; OS* [bW] S _ P ; 

03«=[» q J] s _ p [bJ q <] s _ p ; 

OLjli ee ^{fiV-ltflv-AtFrtv-A + F?iv-A[VF-7 ? >_ A } ; (2.1) 

with q being a strange or a down quark, and i, j colour indices. The last operator, as well as similar 
1/M corrections OSjl q and 03jl q for the OS q and 03 q operators, are required at 6 (Aqcd/M). 

The continuum matrix elements {OX)(jj,) MS = (B^|OX|B^) MS (^) of the operators OX = OL q , 
OS q , 03 q entering in the SM formulae, are related to those evaluated via lattice simulations by 
a perturbative one-loop matching relation through &{a s ), & (Aqcd/M) and &(a s /{aM)). The 
matching relations mix, already in the continuum, the four-fermion operators in ( |Q| ) -see [Q, ||] for 
the explicit expressions. 
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Table 1: Simulation parameters for the coarse (first four sea masses) and fine lattices (last line) for B S Q 

(m v q aL /mf ys - = 1) and B d Q (m v q aL = m s d ea ). 



The bare hadronic matrix elements are obtained by numerically evaluating the three-point and 
two-point correlation functions 

cW{tuti) = £ <0|1> 5? (*i,/i) [Q] (0)^t(^ 2 ,-/ 2 )|0) 

C( B )(0=£(0|<I> 59 (f,0<I>t(0,0)|0) (2.2) 



with <J> B - ? (x,?) = b{x,t)y=,q{x,t) and O any of the four-fermion operators in (1A). The simula- 
tions are performed on MILC configurations with Nj ■ = 2 + 1 sea quarks. The valence b fields 
are described by the NRQCD action improved through 0(1 /M 2 ), &{a ) and leading relativistic 
G{\ /M 3 ) [Q], while the light valence (and sea) quarks are staggered asqtad fields [Bt]. An improved 
gluon action is also used to further reduce discretization errors. 

The action parameters are fixed via light and heavy-heavy simulations, in particular the valence 
b and s quark masses are tuned to give the physical values of the T and K mesons. The different 
parameters in the simulations are collected in Table [lj 



2.1 Mixing parameter for Bq mixing 

Our work published in [||] analyzes the B® mixing parameters for two ensembles of MILC con- 
figurations with (m^" = mf Ja )/m s = 0.25,0.50 and a = O.llfm (coarse lattice). This corresponds 
to the first two entries in Table |l| with m v q al - /m% hys ' = 1. 

The results obtained for the mass and width differences when using these parameters in the 
SM expressions are 

AM i = 20.3(3.0)(0.8) j p^ 1 and Ar. v = 0.10(3)p^ 1 , (2.3) 

which agree with experimental results within errors. The first error in AM S , which is the dominant 
one, is from the lattice determination of fsfiB s through the definition 

(OL)jg = (B s \OL\B s )™ = \f 2 B B Bs {n)Ml , (2-4) 

and the second one is an estimate of the error from |V r *V^| and rn t . This 15% lattice error is dom- 
inated by a 9% statistics+fitting error and a 9% uncertainty associated with higher order operator 
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matching. The large statistical errors are due to the fact that the simultaneous fits of two-point 
and three-point functions are unstable and we need to constrain the two-point parameters using the 
values obtained in fits to only two-point correlators. 

The stability of the fits can be improved by using smearing techniques that reduce the overlap 
with excited states. We checked that the statistical+fitting error can be reduced from 9% down to 
as low as 2% by smearing the heavy quark in the two-point functions and further improvement is 
achieved by smearing also in the three-point functions, as described in the next section. 



3. New results for B s and Bq mixing parameters 

We have generated two-point functions with both local and smeared sources and sinks, using 
a smearing we call IS since it takes an exponential form. Our three-point functions are local at the 
source (the site of the 4-quark operator), and have both local and smeared sinks at either end, with 
the same smearing as in the two-point case. This reduces the statistical+fitting error in our analysis. 



The general definition of these two-point and three-point functions is given in (2.2). 



In addition to the matrix elements relevant in the determination of Bq mixing parameters, we 
have also calculated those corresponding to Bf } mixing in full QCD. The s and b valence quarks 
masses are the physical ones, while the d valence quark mass is the same as m s d ea for any ensemble. 
Two different lattice spacings have been studied, the MILC coarse lattice (a = 0.12) and the MILC 
fine lattice (a=0.09). On the first one we have the correlation functions calculated for four different 
values of the light sea quark masses and on the second one, so far we have results only for one 
light sea quark mass. The parameters of the simulations, quark masses, number of configurations, 
number of time sources, etc, are shown in Table |l[ 

We have not analyzed yet the l/M corrections for all the data collected in Table 0, so the 
results presented in these proceedings are only coming from the dominant contribution in the 1 /M 
expansion. We are also still working on the fits with the lightest sea mass on the coarse ensemble, 
m q = 0.005, and results for this point will be presented elsewhere [|]]. 

3.1 Reduction of statistical+fitting errors 

The use of several time sources and smearing greatly reduce the statistical errors as can be seen 
in Figure |I[ In that Figure, as an example of that reduction, we compare the results in our previous 

paper ^ for fs s \l Bs s {GeV) with our new results incorporating smeared correlation functions in 

the fits and new data. fg s \J BB s (GeV) is plotted as a function of the light sea quark mass over the 
physical strange quark mass, m q /m s , and the errors are only statistical. 

With these new data we are able to get stable simultaneous fits for two-point and three-point 
correlation functions without any constraints in the two-point function parameters. With stable we 
mean that the central values, errors and % 2 jndof do not change when we add more excited states 
in the f uncti onal forms to be fitted. The re sult is a reduction of statistical errors from 4.5% to 1-2% 
in f Bs ^B\ s {GeV), and similarly for f Bd ^B\ d . 

Another technique that could reduce further the size of statistical errors is the use of random 
wall sources for the light propagators. We have already checked that the statistical errors in the B*' 
two-point parameters are improved by a factor of two, comparing results for the same heavy-light 
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Figure 1: Bs s {GeV). Errors are only statistical. 



correlators we are using here but with HISQ [fT]] (Highly Improved Staggered Quarks) light valence 
quarks, with and without random wall sources -see ^ for more details about using random wall 
sources in heavy(NRQCD)-light(HISQ) correlators. Further study is needed to find how the use of 
this kind of source affect the three-point function parameters relevant for B° mixing. 

3.2 Calculation of the ratio ^ 

Some of the errors affecting the calculation of fy q \J^B q will cancel almost completely and 
others partially in the ratio t, = ^ S ^^L . In Figure || we show values for this ratio multiplied by 



the square root of the masses of the Bq and Bq mesons, 



(3.1) 



together with the ratio O s /<& ? 



plotted as a function of m q /m^ hyi '' , where m q = «j 

The errors for X s /X q in Figure (^), which are only statistical, are larger than those for <& s /<& q 
because we have not yet taken into account the correlations between the data in the numerator and 
denominator in this ratio. We expect to reduce this error to less than 2% when these correlations are 
included (the current plotted values have 2.5% errors). Another error that should be significantly 
reduced is that for the fine lattice point since we have not yet included all of our data. 



fn s jM Bs 



>M R , 



, without the bag parameters from \w. The results are 



.valence 



m„ 
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as a function of the light valence quark mass in the 



denominator. Errors are only statistical in all the quantities plotted. 



The statistical errors are not the only ones to be reduced by taking the ratio. Discretization, 
relativistic and higher order operator matching will affect fs s \JBb, and fg d in the same way 
and largely will cancel in the ratio. One expects their effects to come in at the level of the cor- 
responding error in fg q \JBb {I times a{m s — mj) or (m s — md) /Aqcd- The results for /bJ fB d are 
nearly unchanged when adding one-loop and l/M corrections [^] and we expect something similar 
here. We have already checked that the difference between tree level and one-loop results is less 
than 1%. The scale a~ 3 uncertainties, that lead to a 5% error in /g Bs q , do not affect the ratio ^ . 

The next step in our calculation will be to carry out a chiral extrapolation of these results to the 
physical point including the effect of taste-changing errors, to account for the remaining systematic 
in the calculation and remove the dominant light discretization errors. 



4. Summary and future work 

We have calculated the mixing parameters in the B® and B° d systems for two different lattice 
spacings and five different light quark masses. The statistical errors have been reduced from our 
previous work by a factor of 2-3, so statistics is no longer a dominant source of uncertainty in 
the calculation of f^Bs q - The largest error is now the uncertainty associated with the perturbative 
matching, that is also reduced from 9% to 6.5% by simulating on finer lattices. Further reduction 
of this source of error, as well as discretization errors, would also be possible by the use of MILC 
superfine lattices. 
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We also give preliminary results for the ratio E, versus m q /m s , where many theoretical uncer- 
tainties are partially or completely cancelled between denominator and numerator. 

The analysis of l/M corrections and results for mdjm s = 0.125 and at least one other light 
quark mass on the fine lattice, will be presented in a forthcoming publication We are also 
exploring different smearings and better fitting approaches to further reduce the statistical errors. 
In particular, we are getting promising preliminary results using random wall sources for the light 
propagators. 

Once other sources of errors have been reduced, we need to perform a chiral extrapolation 
of the /bVBb and £ results incorporating light discretization uncertainties (taste-changing errors) 
and perturbative errors. We will also be able to perform a continuum extrapolation, since we have 
results for two different values of the lattice spacing. 

Other talks on unquenched calculations of Bo mixing parameters in this conference can be 
found in M\. 
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